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The distribution of spacings between the lasing frequencies for an ensemble of random lasers in 
the two-mode regime was computed. The random lasers are implemented as open chaotic cavities 
filled with an active medium. The spectral properties of the passive cavities are modeled with non- 
Hermitian random matrices. The spacing distribution is found to depend on the relation between 
the gain-profile width and the mean spacing of the passive-cavity modes. The distribution displays 
mode repulsion and, under certain conditions, agrees with the Wigner surmise. The role of mode 
competition is discussed. 

PACS numbers: 42.55.Zz, 05.45.Mt 

I. INTRODUCTION II. THEORY 



The term "random laser" usually refers to lasing sys- 
tems based on disordered materials or substantially open 
wave-chaotic resonators [1]. They are distinguished from 
conventional lasers by leaky modes (because of the ab- 
sence or openness of the resonator) and almost random 
field distributions (due to disorder or chaotic shape). 
Coherent lasing in disordered materials has been ob- 
served via emission spectra [2] and photon-count mea- 
surements [3]. The role of strong localization of light 
in disordered lasers remains an open theoretical ques- 
tion [4, 5]. In the recent work on the chaotic-laser the- 
ory [6, 7], the standard laser models [8, 9] were ex- 
tended to cavities with spectrally overlapping modes and 
equipped with the ideas of quantum chaos [10]. The spa- 
tial structure of lasing modes was also considered [11]. 

A substantial aspect of the probabilistic description 
of random lasers is the mode statistics. For example, 
the mean number of lasing modes in weakly [12] and 
strongly [13] open resonators, as well as the variance [14], 
were calculated with the help of random-matrix the- 
ory [15]. The present study was partly motivated by 
recent experiments in a porous-GaP laser [16]. The ob- 
served distribution of spectral mode spacings could be 
well fitted with the Wigner surmise for the Gaussian or- 
thogonal ensemble (GOE) of random matrices [10]. As a 
first step toward understanding of the multimode regime, 
here I propose a theory of spacing statistics for two-mode 
lasing. The calculation is based on the random-matrix 
approach of Refs. [13, 14]. The results point to the 
spectral repulsion of lasing modes on the scale inherited 
from the passive cavity. A different repulsion mechanism, 
requiring inhomogeneous broadening, was suggested in 
Ref. [17]. 
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A. Laser equations 

As a model of a random laser we adopt an open res- 
onator of irregular shape, having almost random eigen- 
functions [18]. The lasing takes place due to a pumped 
active medium inside the resonator. Both the field and 
the active atoms interact with the environment. The 
environment variables can be eliminated from the equa- 
tions of motion by inclusion of the effective damping and 
noise terms [8]. The number of photons in a mode A; is 
controlled via the bosonic operators a k and a\. The oc- 
cupation of the two active levels of an atom p will be de- 
scribed by the pseudospin-1/2 operators a p and s p = sj 
satisfying the commutation relation [<7 p ,s p ] = a p . The 
operator a p transfers the electron from the upper to the 
lower level and s p yields the population difference be- 
tween the levels. We will work in the classical approxi- 
mation, replacing operators with c numbers and neglect- 
ing the noise. Then I k = \a k \ 2 will be proportional to the 
mode intensity and a p will characterize the atom's polar- 
ization. The coupled classical equations of motion (equiv- 
alent to the Maxwell-Bloch equations) take the form 

a k = - (iu k + Kk) a k + g^(f>* k (ry)oy, (1) 

p' 

& p = - (iv + 7x) gp + 2g <f> k > (r p )a k >s p , (2) 

k' 

s P = 7|| ( s ~ s p) -9^2 \$ k ' ( r p) a fe ,cr p +c.c] . (3) 

k' 

Here u) k — in k are the complex frequencies of passive 
modes in the cavity, v is the atomic transition frequency, 
-f± and 7|| are the polarization and inversion decay rates, 
and S specifies the pump strength. The atom-field cou- 
pling is measured by the parameter g ~ dy^rns/h, where 
d is the dipole moment for the atomic transition. All 
modes k are assumed to be linearly polarized in the same 
direction. The field amplitude, evaluated at the atom 
position Tp, is described by the normalized wavefunc- 
tion (f>k (r) . As an eigenfunction of a non-Hermitian oper- 
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ator, it has an adjoint eigcnfunction 4> k (r) (Cf. Ref. [19]). 
The functions (f>k and <j>k' are orthogonal for k ^ k' 
(biorthogonality) . Equations (l)-(3) are written in the 
rotating-wave approximation. It neglects antiresonant 
products of the type a,k<J p , oscillating with a double op- 
tical frequency. 



B. Lasing frequencies 



To solve Eqs. (l)-(3), we assume that each las- 
ing mode dfc(t) = afc(i) exp(— ififci) oscillates with a 
constant frequency fife and a slowly varying ampli- 
tude afc(t). Following a standard procedure [9] we 
expand er p (t) = ^ fc a p k(t) exp(— ififcf) and s p (i) — 

= J2kk> i s pkk' (t) exp [-i (fife - fife/) i] + c.c.} and ne- 
glect all other possible oscillating terms. Furthermore, 
only contributions up to the order of a| in a p and a k in 
s p are kept . Under the assumption | ctk/ctk | -C 7j_ , 7|| , the 
atomic variables can be eliminated from the equations of 
motion. The resulting equations for afe, 



If only one mode is lasing, its frequency is [9] 



fi 



(1) _ U) k + UKk/jA 



1 + K k /j± 



(8) 



The positive-intensity condition requires that the linear 
gain for this mode exceeds the linear absorption. Hence, 
the actual first lasing mode, labelled k = 1, without loss 
of generality, is the one with the lowest threshold 
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In other words, the first lasing mode k = 1 emerges at 
the pump level G = G±. Unlike the frequency of 

the second lasing mode depends on the mode intensities. 
In the results for the spacing statistics below, this fre- 
quency is always taken at the threshold. If k ^ 1 is the 
second lasing mode, Eq. (4) provides a pair of complex 
equations, with an unknown I\ and I k — > + . They yield 

(2) 

a single nonlinear equation for the frequency tt k at the 
threshold, 



"! = <•:■ | Pk - qklk - ^ r kk'h> 

k'^k 



(4) 



contain linear and nonlinear (intensity-dependent) terms. 
The coefficients arc 



p k = i (Ofe - Lo k ) - Kfe + GV k , 
qk = 2cBk kkk T>kCk-, 



r k k> = cD k 



^Bkkk'k'^k' + Bkk'kk'T^kk' (Dk> + ^fc) 



(5) 
(6) 



(7) 



where B Mmn ~ Vj dr <p* k (r) <pi(r) <p m (r) ^*(r), V k = 

[l-i(fi fc -i/)/7±r\ V\ k , = [l-i(n fc -fi fc ,)/7||] _1 , 
C k = RePfe, G = 2g 2 SM/ lA _V, c = Ag^SM/^l 7||^ 2 , 
J\f is the number of atoms, and V is the cavity volume. 
The linear absorption Kk is counteracted by the linear 
gain GCk , which has a Lorentzian profile of halfwidth 7^ . 
The chaotic wavefunctions for different modes behave like 
independent Gaussian random functions, on scales larger 
than the wavelength [18]. This allows one to approx- 
imate the correlation parameters Bkkk'k' — Bkk'kk' — 
1 + 2J fcfc / [12, 13]. 

In a steady state a k = 0, the bracketed portion of 
Eq. (4) must vanish for all k, such that a.k 7^ 0. This 
makes a system of complex algebraic equations, from 
which the fife's and Ik's of the lasing modes can be deter- 
mined. To use this procedure, however, it must be known 
a priori which modes have nonzero intensity. If this is 
not the case, one has to assume that certain modes are 
lasing, solve the system of equations, and then verify that 
all found intensities are positive and the solution is sta- 
ble. In view of the complexity of the problem, I restricted 
the present study to the case of two-mode lasing. 
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The stability of the two-mode solution can be checked 
as described in Ref. [8]. Of all stable candidates, the 
genuine second lasing mode, labeled k = 2, minimizes 
the threshold, i.e., 



G 2 = min (d+QC^Rk) ■ 

stable fc^l V / 



(12) 



The first-mode intensity at the second-mode threshold is 
Ii = R2/C. In order to estimate the importance of the 
nonlinear effects (mode competition), we can also define 
the would-be second mode neglecting the I\ dependence 
in Eq. (4). This mode, labeled k = 2', has the threshold, 

simply, G2' = minfe^i (^Kk/C^^j and the frequency fi^. 



C. Random-matrix model 

The working equations (8)-(12) use the eigenfrequen- 
cies Wfc — ink of the passive cavity as an input. In or- 
der to describe statistical characteristics of an ensemble 
of (nonidentical) chaotic lasers, the passive spectra can 
be modeled with the eigenvalues of non-Hcrmitian ran- 
dom matrices. Henceforth we will formally set v = 0, 
i.e., the frequencies will be measured with respect to the 
atomic frequency. Following Rcfs. [13-15], we associate 
with each cavity an L x L matrix u> — ij. Here, the real 
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symmetric matrix Co is chosen from the GOE and 7 is a 
fixed diagonal matrix with M <C L positive and L — M 
zero eigenvalues. The eigenvalue de nsity of Co obeys the 
Wigncr semicircle law p(uo) = n^ 1 yl — w 2 /4 (in dimen- 
sionless units) , which remains approximately valid for the 
real parts of the eigenvalues of Co — 17. The integer M is 
interpreted as the number of spectral bands of the out- 
side field, to which the cavity field is coupled [7]. The 
bands, called the coupling channels, are a consequence of 
the wave quantization due to the finite size of the cav- 
ity opening. They are similar to conductance channels 
of a quantum dot connected to ballistic leads. We will 
use the model of equivalent channels, when all nonzero 
eigenvalues of 7 are equal to some 7 > 0. According to 
random-matrix theory, the coupling strength is charac- 
terized by the parameter 2np (lo) / (7 + 7~ 1 ) [15]. Hence, 
7 between and 1 spans the whole coupling range from 
weak to strong. In order to exclude the effects of vari- 
able density p(oo), only the eigenvalues near the top of 
the Wigner semicircle were taken into account in the nu- 
merical simulations. That is, for each random matrix 
Cj — i-y, I used L w 0.36L eigenvalues with the smallest 
thereby allowing for a 4% variation of p(wfe). 



III. NUMERICAL RESULTS AND DISCUSSION 

Numerical results for the mode-spacing statistics are 
shown in Figs. 1 and 2. The first two lasing frequen- 
cies O^ and O^, as well as the second frequency O^ 
without mode competition, were computed for each ran- 
dom matrix in an ensemble. The plots show the spacing 
distributions P(AO), P' (AO') (the primed function is 



not a derivative), and p (Aw), where AO = 
AO' = 



0« - 2 



(2) 



0« 



, while Ato are the nearest-neighbor 

spacings for the real parts LOk of the passive eigenfrequen- 
cies. (The spacings in the figures are scaled with their 
respective mean values.) p (Aw) for 7 <C 1 is close to the 
Wigner surmise 



p w (Auo) m — Alo cxp ( — — A. 



(13) 



which displays eigenvalue repulsion in the GOE. For 
stronger coupling, the eigenvalues get spread in the com- 
plex plane, and the crossing probability for the real parts 
increases. According to the numerical data (Table I), the 
average spacings AO and AO' are of order of the gain- 
profile halfwidth j±. We will consider the three regimes 
of j± being much greater, of order, or much smaller than 
the passive mean spacing Alo « 1/Lp(0). The inversion 
decay rate 7|| is taken equal to j±. 

1. 71 » Aw [Figs. 1(a) and 2(a)]. Although the dis- 
tributions P (AO) and P' (AO') are quite similar, the 
mode competition still influences the lasing frequencies. 
In principle, this influence is twofold: (1) the indices 
k = 2 and k = 2' refer, in general, to two different pas- 
sive modes and, (2) even if the modes are the same, the 
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FIG. 1: Distributions of spacings between the first two lasing 
modes for an ensemble of 6,000 non-Hermitian random ma- 
trices of size L = 200 with M = 6 coupling channels and cou- 
pling parameter 7 = 0.1. In each matrix, L = 72 eigenvalues 
closest to the top of the Wigner semicircle were taken into ac- 
count. The average damping TcjT = 3.0 x 10~ 3 was determined 
numerically. Shown are the distributions P (AQ) (bold solid 
line), P' (AQ') (circles) with the mode competition neglected, 
p (Alo) (pluses) for the passive frequencies, the Wigner sur- 
mise pw (Alo) (thin solid line), and the approximation (16) 
for P' (AQ') (dashed line). The spacings are computed in 
units of the mean values for each distribution (Table I). The 
gain-profile halfwidth -y± is 10 _1 (a), 10^ 2 (b), and 10~ 3 (c). 
The inversion decay rate j\\ = j±. 
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AQ 


AQ' 


Alo 




10~ 3 


0.0049 


0.0043 




0.1 


10~ 2 


0.017 


0.015 


0.016 




10- 1 


0.102 


0.079 






10- 3 


0.0020 


0.0016 




1.0 


10~ 2 


0.017 


0.014 


0.015 




10- 1 


0.11 


0.087 





TABLE I: Mean frequency spacings AQ between the lasing 
modes (AQ' neglecting the mode competition) and Alo be- 
tween the passive-cavity modes for different 7 and 7x = 7y . 



is different from O^ 1 ' [Eq. (10)]. It can 



frequency O^ 2 "* 
be deduced from Eq. (10) and was also checked numer- 
ically (not shown) that, in the present case, the second 
factor is insignificant. On the scale AO > AO, the pas- 
sive distribution p (Alo) is irrelevant for P (AO), because 
the passive modes corresponding to k = 1, 2 are not the 
nearest neighbors. The passive modes with close Wfe's are 
less likely to become lasing modes, since they may lie far 



4 



a 
° i 

j5 0.8 



£ 0.6- 



1 

0.8 



0.2 




1 1 1 1 

/ • /* \ 


i 

(c) : 






/ %^ 


- 


— 1 — 1 — — - 


(b) - 




- 


/ 




1 1 

• 

_ • 
• /a 

A* a ~\ 


(a) - 


•// »»V\ 

7/ 

/ , i , i '** 





12 3 4 

mode spacing [in mean spacings] 

FIG. 2: Same as in Fig. 1, but with 7 = 1.0, 5 = 4.2x 10" 



from each other in the complex plane, i.e., one of them 
can be substantially damped. This explains the falloff 
of P (Ail) as Ail — ► 0. For small 7, the repulsion of bas- 
ing frequencies is guaranteed by the Wigner repulsion of 
the passive modes. 

An approximate expression at Ail' 3> Alu for the dis- 
tribution P' (Afi') without mode competition can be de- 
rived from a simple probabilistic model as follows. Of all 
passive eigenvalues u)k — i^k, the two lasing modes have 
the smallest Kk = Kk (l + w^/7^) [here we can approxi- 
mate il^ « uik, since only small n k are relevant]. Let us 
divide the lu axis into 3> 1 intervals, each containing 
L/Nu 3> 1 eigenvalues. These conditions ensure that the 
two lasing modes arc unlikely to belong to the same in- 
terval and that the Kk distributions in different intervals 
are uncorrclated. The smallest Kk within an interval is 
distributed according to P gap (k) = (dn/dn) e~ n ^ [as- 
suming that the number of eigenvalues between k and 
k + dn obeys the Poisson distribution with the aver- 
age dn (k)]. Here n (k), the mean number of eigenvalues 
with Kk < k, has a small-/* asymptotics k~ [15] (valid 
for arbitrary 7). The distribution of the smallest Kk is 
P^ p (k) = P gap (k) dn/dn, where is the central fre- 
quency of the interval. The two lasing modes lie in the 
intervals i and j with the probability 



J dKdK'P£ p (K)P£ p (K') 
/>max(/5,/5') 

Jo 



dK"P£ p (K") 



■ (14) 



Using an argument similar to that of Ref. [10] (Sec. 5.6), 
the product in the second line can be estimated as 



exp j-0 Lmax(K,K') M/2 j. Thus, only k < L~ 2 / m 
is relevant. The spacing distribution is 

P'(Ail') oc J dujP£" +AQ ' oc J dto (£, 



r \ M / 2 



(15) 

where C u = (l + w 2 /7^) and the integration was ex- 
tended to infinity. For M = 6 (see the figures), 

o,/aom 32 a; 4 + 24.x 2 + 336 Ail' . . 

P'(Ail')^ — , xee— . (16) 



3^7± (x 2 + 4)' : 



7-L 



2. 7j_ < Aw [Figs. l(b,c) and2(b,c)]. The lasing modes 
are likely to emerge from the passive nearest neighbors. 
Thus, the spacing distribution approaches the Wigner 
surmise (13). At small the mode competition is man- 
ifested mostly via the frequency shift (10), while k = 2 
and k = 2' is usually the same mode. This frequency 
shift is responsible for the vanishing P (Ail) as Ail — > 0, 
even when P' (Ail') does not vanish [Figs. 1(c) and 2(c)]. 

In other words, when il^ S e ^ s close to il^\ il^ i s rc ~ 
polled from these two frequencies. It is expected that 
the distribution P 1 (AO') without the mode competition 
is close to p(Alu). In particular, when -f± <C k& ~ Auj 
[Fig. 2(c)], the frequencies (8) can be approximated as 



il 



(i) 



ujk"fA_/Kk- This makes C 



(i) 



1 + (LUk/Kk) 

i.e., the gain profile becomes independent of 7j_. When, 
on the other hand, 7_l ^ Kfe <C Aw [Fig. 1(c)], the two 
modes with the smallest \ilk\, not the smallest Kk, are 
lasing. Since k^s for these modes can be substantially 
different, Ail' is not proportional to Auo [Eq. (8)]. This 
explains the difference between P' (Ail') and p (Au) in 
Fig. 1(c). 



IV. CONCLUSION 

The spacing distribution for the first two lasing modes 
depends on the relation between the gain-profile width 
and the mean spacing of the passive-cavity modes. For a 
very wide gain profile, the passive-spacing statistics is ir- 
relevant, except for very small spacings. When these pa- 
rameters are of the same order, the spacing distribution 
is approximated by the Wigner surmise. The nonlinear 
frequency shift due to the mode competition provides for 
the mode repulsion, even when the unperturbed frequen- 
cies may approach each other. 
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